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A model for the asteroid resurfacing by regolith convection is built to estimate its timescale. In the model, regolith convection is 
driven by the impact-induced global seismic shaking. The model consists of three steps: (i) intermittent impact of meteoroids, (ii) 
impact-induced global vibration (seismic shaking), and (iii) vibration-induced regolith convection. In order to assess the feasibility 
of the resurfacing process driven by regolith convection, we estimate the resurfacing timescale as a function of the size of a target 
asteroid. In the estimate, a set of parameter values is assumed on the basis of previous works. However, some of them (e.g., seismic 
quality factor Q, seismic efficiency rj, and seismic frequency /) are very uncertain. Although these parameter values might depend 
on asteroid size, we employ the standard values to estimate the representative behavior. To clarify the parameter dependences, we 
develop an approximated scaling form for the resurfacing timescale. According to the estimated result, we find that the regolith- 
convection-based resurfacing timescale is shorter than the mean collisional lifetime in most of the parameter uncertainty ranges. 
These parameter ranges are within those reported by previous works for small asteroids. This means that the regolith convection 
can be a possible mechanism for the asteroid resurfacing process. 
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1. Introduction 

The history of surface evolution of small asteroids is 
recorded on their surface terrain. Impacts of various-sized me¬ 
teoroids are the most important processes to discuss the asteroid 
surface evolution. For instance, regolith grains are generated 
and mixed by repeated impact events. This gardening effect 
keeps the surface fresh and active (e.g., Melosh (2011)). 

However, recent planetary explorations have revealed a rich 
variety of surface terrains that cannot be understood only by 
the impact cratering and resultant gardening effect. For in¬ 
stance, asteroid Eros, which was explored by NASA’s space 
probe NEAR Shoemaker from February 2000, showed the in¬ 
teresting characteristics. Based on the observational result, As- 
phaug et al. (2001) reported the surface blocks distribution that 
cannot be explained by the potential source craters. To explain 
the distribution, they discussed the possibility of the size sort¬ 
ing in the vibrated regolith bed. More recently, Miyamoto et al. 
(2007) and Tancredi et al. (2015) also reported that the sur¬ 
face regolith of asteroid Itokawa, which was explored by the 
Japanese space probe Hayabusa, could migrate and be segre¬ 
gated. In addition, impact craters on the surface of Itokawa are 
very vague (Saito et al., 2006; Hirata et al., 2009). The analysis 
of tiny samples returned from Itokawa also supports regolith 
migration. Tsuchiyama et al. (2011) observed some round- 
shape particles that might be worn away by friction between 
regolith grains. Moreover, the cosmic-ray exposure (CRE) age 
of Itokawa’s returned particles are measured by Nagao et al. 
(2011) and Meier et al. (2014). According to these studies, the 


upper limit of CRE age is approximately 8 Myr and its proba¬ 
ble value is 1.5 Myr. This value is relatively young compared 
to the CRE age of most of the EL chondrites. All the above- 
mentioned observations and findings suggest that the regolith 
migration process plays a crucial role in the surface evolution 
on small asteroids such as Eros and Itokawa. 

Eor the driving force of regolith migration on 
Itokawa, Miyamoto et al. (2007) followed the idea of As- 
phaug et al. (2001); they considered the size segregation of 
vibrated regolith layer. They supposed that the impact-induced 
global seismic shaking could result in the convection of re¬ 
golith grains. Moreover, the size segregation of regolith grains 
could also be driven by regolith convection. Although the 
impact-induced seismic shaking has not been directly observed 
on asteroids, Richardson et al. (2004) and Richardson Jr 
et al. (2005) developed a model of the global seismic shaking 
on asteroid Eros to investigate the crater erasure process. 
Recently, Garcia et al. (2015) has also computed seismic wave 
propagation on asteroids in kilometric and sub-kilometric size 
ranges. Using the global seismic shaking model of Richardson 
et al. (2004) and Richardson Jr et al. (2005), Miyamoto et al. 
(2007) showed that the regolith convection induced by global 
seismic shaking can be caused even by small-scale impacts. 
Eor instance, 1 cm impactor with the impact velocity Vi = 4 
km s“* is sufficient to achieve the global seismic shaking on 
Itokawa because the Itokawa’s surface gravity is extremely low, 
lO^'* m s“^ (Abe et al., 2006; Eujiwara et al., 2006). Thus, the 
effect of collective granular motion on the surface of asteroid 
could be very important for properly considering the surface 
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evolution (Asphaug, 2007; Katsuragi, 2016). 

In fact, granular convection can readily be observed when 
the granular matter is subjected to the vertical vibration in the 
laboratory experiment (e.g., Yamada and Katsuragi (2014)). In 
granular convection induced by vertical vibration, the most im¬ 
portant parameter is a ratio between the maximum vibrational 
acceleration and gravitational acceleration, F = Ao(27r/)^/g, 
where, Aq, / and g are the amplitude of vibration, its fre¬ 
quency, and gravitational acceleration, respectively. The on¬ 
set criterion of granular convection can be approximately writ¬ 
ten as r > 1 (e.g., Garcimartm et al. (2002)). In addition, if 
the vibrated granular matter consists of polydisperse grains, the 
vertical size segregation of grains can also be observed. This 
vibration-induced size segregation is called Brazil nut effect 
(BNE). Granular convection is considered a plausible reason for 
the BNE in the vertically vibrated granular bed (Knight et al., 
1993). 

Recently, the gravity dependence of granular convective mo¬ 
tion has been investigated on the basis of laboratory exper¬ 
iments (Murdoch et al., 2013b,c,a; Guttler et al., 2013; Ya¬ 
mada and Katsuragi, 2014) and numerical simulations (Tan- 
credi et al., 2012; Matsumura et al., 2014). Murdoch et al. 
(2013b,c,a) investigated the gravity dependence of the veloc¬ 
ity of shear-driven convective motion by using the parabolic 
flight method. They reported that the convective velocity is sig¬ 
nificantly reduced under the microgravity condition while it is 
enhanced under the high-gravity condition. Guttler et al. (2013) 
also performed an experiment of convection-driven BNE using 
the vibrated granular bed on parabolic flights. They found that 
the rise velocity of an intruder in the vibrated granular bed is 
roughly proportional to the gravitational acceleration g. Tan- 
credi et al. (2012) numerically computed the threshold velocity 
Vth to induce BNE and reported Vth ~ It should be noticed 
that, however, v* is the threshold velocity above which segre¬ 
gation starts. Matsumura et al. (2014) numerically investigated 
the convection-driven BNE and found that the rise velocity of 
an intruder in the vibrated granular bed is proportional to g^'^. 
In order to evaluate the gravity dependence of granular convec¬ 
tive velocity, the scaling method has been used in Yamada and 
Katsuragi (2014). The researchers investigated the convective 
velocity under the steady vertical vibration on the ground. Us¬ 
ing an experimentally obtained scaling relation, they found that 
the convective velocity is scaled by While the intruder’s 
rise velocity (BNE rate) and convective velocity are different 
quantities, Yamada and Katsuragi (2014) is roughly consistent 
with Guttler et al. (2013). However, the reason of difference be¬ 
tween Guttler et al. (2013) and Matsumura et al. (2014) remains 
unclear. In any case, all these experimental and numerical re¬ 
sults suggest that the granular convection can occur with very 
low velocity even under the microgravity condition. 

In general, the characteristic timescale for gravity-related 
phenomena is elongated under the microgravity condition. Eor 
instance, the penetration timescale of an intruder into a granu¬ 
lar target increases under the low-gravity condition (Altshuler 
et al., 2014). Eor granular (regolith) convection, the character¬ 
istic timescale would become very long under the microgravity 
condition because the convective velocity, which is inversely 


proportional to the timescale, is significantly reduced under the 
microgravity condition as mentioned above. Therefore, it is not 
obvious whether the low convective velocity can really resur¬ 
face the various-sized S-type asteroids covered with regolith 
within the reasonable timescale. The estimate of timescale is 
necessary to assess the feasibility of the resurfacing by regolith 
convection. The main purpose of this study is to build a model 
of resurfacing process induced by regolith convection for aster¬ 
oids covered with regolith. Using the model, timescale of the 
resurfacing can be estimated. Then, the resurfacing timescale 
should be compared with other timescales such as the mean col- 
lisional lifetime. Based on the comparison, we discuss the fea¬ 
sibility of the resurfacing process driven by regolith convection 
on the surface of asteroids. 

This paper is composed of following sections. In the next 
section, we build a model of the resurfacing process. Section 3 
presents the estimate of resurfacing timescale. In addition, the 
scaling analysis for the estimated timescale is discussed in this 
section. In section 4, we compare the resurfacing timescale 
with other timescales and verify the feasibility of the regolith- 
convection-drivenresurfacing. Section 5 mentions conclusions. 


2. Model 

2.7. Outline of the resurfacing model 

In order to estimate the resurfacing timescale, we make a 
simple model of the asteroid resurfacing process driven by re¬ 
golith convection. The regolith convection is caused by impact- 
induced global seismic shaking. In the model, we divide the 
resurfacing process into three steps as follows: 

(i) Impact: Impactors intermittently collide with a target as¬ 
teroid, 

(ii) Vibration: Each collision leads to a seismic shaking, 

(iii) Convection: The seismic shaking induces the regolith con¬ 
vection on the asteroid. 

Eig. 1 shows the schematic illustration of the model. Repeat¬ 
ing the series of three steps, the asteroids can be resurfaced by 
regolith convection. 

In the convection step, we assume that the iteration of 
small convective migrations constitutes the well-ordered (re¬ 
producible) convective roll structure. In Eig. 1, each red ar¬ 
row corresponds to the migration length driven by each impact 
event, and they collectively construct the well-ordered identical 
convective roll. Although one may feel that this assumption is 
not so plausible, some experimental works have reported that 
the intermittently tapped granular bed can construct the well- 
ordered convective motion or BNE (Philippe and Bideau, 2003; 
Ribiere et al., 2005; Tancredi et al., 2012; likawa et al., 2015). 

The details of steps (i), (ii), and (iii) are written in sec¬ 
tions 2.1.1, 2.1.2, and 2.1.3, respectively. Einally, we unify 
these three steps and compute the resurfacing timescale. The 
unification is described in section 2.2. 
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(i) Impact Surface migration length 



Figure 1: Schematic illustration of the model for the asteroid resurfacing in¬ 
duced by regolith convection. Intermittent regolith convection results in the 
resurfacing of asteroids through three steps: (i) impact, (ii) vibration, and (hi) 
convection. Each red an'ow’s length conceptually shows the convective mi¬ 
gration length per impact event, /. We assume that the convective roll size A 
corresponds to the surface migration length of the convecting regolith. 

2.1.1. Impact frequency 

The frequency of impacts per year A^p is a necessary quantity 
to build the model. A^p is estimated from the population of small 
objects A^i of diameter D^. To estimate Ni for the Main Belt 
Asteroid (MBA) and Near-Earth Asteroid (NEA), there have 
been two models based on the dynamical evolution of the main 
belt population: O’Brien and Greenberg (2005) (called here¬ 
after OBG population) and Bottke Jr et al. (2005) (called here¬ 
after BAL population). The range of estimated impactor size in 
OBG population is broader than that of BAL population, and 
both are in good agreement with the currently observed popu¬ 
lation in the main belt. Only the minor differences can be rec¬ 
ognized between these two models. Here, we employ the OBG 
population to compute Np. Eig. 2 shows the cumulative number 
distribution of potential impactors Ni,cum(^ D') for MBAs and 
NEAs, estimated by OBG population. 

By numerically differentiating M,cum, we can obtain the in¬ 
cremental distribution of Ni of small bodies with diameters be¬ 
tween D', and D',,,, where k is the index of the binned data 

> D[^. The relation between Ai_cum(S D'n) and Ai(Di ,t) 

is described as 

NfDpk) = M,cum(> D[j) - M,cum(> (1) 

Here, we consider the representative value of size Dik - 
D'^k+i^ik because the OBG population data adapted logarith¬ 
mic bin size for D'^.. 

By multiplying Ni by the annual mean collision probability 
per collisional cross-section Pi and the collisional cross-section 
(A/2+£>a/2)^, we can compute the number of impacts per year 
Np . The specific form of Np is written as 

/ Di Da 

Np(Di, A) = M(A)Ei (Y + y) ’ (2) 



Figure 2: Cumulative number distribution A^i,cum(^ ^[) "vs. impactor diameter 
D'. Red open circles and blue open triangles respectively show the numerical 
data for MBAs and NEAs computed by O’Brien and Greenberg (2005). The 
solid lines are power-law fits. 

where is the diameter of target asteroid. Since this procedure 
is applied to all k data, we omit the subscript k hereafter. In this 
study, we use the constant values Pi - 2.9 x m^^yr^' for 
the collision of MBAs (Bottke and Greenberg, 1993; Bottke Jr 
et al., 1994a) and Pi - 15 x 10“^"^ m^^yr ' for the collision 
of NEAs (Bottke Jr et al., 1994b). Thus, Np can merely be 
defined as a function of Di and (O’Brien and Greenberg, 
2005). Actually, the value of Pi depends on the location and 
orbit of small bodies. Here we neglect such higher-order effects 
and consider the mean value. 

2.1.2. Global seismic shaking 

Next, the strength of seismic shaking E, which corresponds 
to the maximum vibration acceleration normalized to the gravi¬ 
tational acceleration, has to be estimated for each impact event. 
Based on previous reports, we assume that the regolith convec¬ 
tion can be induced by the global seismic shaking satisfying 
r > 1. Moreover, we use the previous experimental result, in 
which the granular convective velocity is scaled by E (Yamada 
and Katsuragi, 2014), to compute the regolith convective veloc¬ 
ity. 

To estimate T caused by meteoroid impacts, we employ the 
global seismic shaking model proposed by Richardson et al. 
(2004) and Richardson Jr et al. (2005). The model is based 
on two assumptions. Eirst, the kinetic energy of the impactor 
is transferred to the seismic energy with the impact seismic 
efficiency factor q. Second, the transferred seismic energy is 
equally distributed to the entire target asteroid by the diffusion 
of seismic energy density. Namely, the homogeneous seismic 
energy distribution is considered as an average picture. The 
distributed seismic energy density attenuates due to the inelas¬ 
ticity of the internal structure of asteroids. The attenuation of 
seismic energy density is characterized by the seismic quality 
factor Q. In addition, spherical shape of the target asteroid is 
assumed for the sake of simplicity. Although the numerical 
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and observational evidences for the validity of this model are 
still very limited, this model offers an analytic way to com¬ 
pute the global seismic shaking on small bodies, and the an¬ 
alytic result was compared with numerical simulations (Garcia 
et al., 2015). From these considerations, F can be calculated 
as (Richardson Jr et al. (2005), Miyamoto et al. (2007) and its 
SOM), 


r(f,A,£>a) 


3/vi 

G 



exp 



(3) 


where pa is the density of the target asteroid, pi is the density of 
impactor, / is the dominant seismic frequency, Vi is the mean 
impact velocity, t is the elapsed time from the moment of im¬ 
pact, and G is the gravitational constant. Fig. 3(a) presents 
an example of F(f) for - 400 m (Itokawa-sized asteroid), 
A = 4 m, and Vi = 5.3 km s ' (Bottke Jr et al., 1994a; O’Brien 
and Greenberg, 2005). Here, the following parameter values 
are also employed; pa = 1900 kg m“^, pi = 2500 kg m“^, 
Q — 2000, rj - 10 and / = 20 Hz (Richardson et al., 2004; 
Richardson Jr et al., 2005; Fujiwara et al., 2006; Miyamoto 
et al., 2007; Yasui et al., 2015). Hereafter, we basically use 
these parameter values as standard ones in the following com¬ 
putations unless otherwise noted. Because we mainly focus on 
the porous rubble-pile or very fractured target asteroids, the rel¬ 
atively low density for the target pa = 1900 kg m“^ is utilized as 
the standard one. For the dominant seismic frequency, Richard¬ 
son et al. (2004) and Richardson Jr et al. (2005) mentioned that 
it falls in / = 10 - 20 Hz for asteroid Eros due to the fractured 
inner structures which act as natural band-pass filter for seismic 
frequency. In addition, they assumed Q = 2000 on the aster¬ 
oid Eros based on the similarity between internal structure of 
fractured S-type asteroids and upper lunar crust. Recently, Gar¬ 
cia et al. (2015) has represented that the dominant frequency 
should increase as the asteroid size decreases. In addition, they 
reported that the relatively high frequency (~ 100 Hz) domi¬ 
nates the vibration on Itokawa-sized asteroid when the mono¬ 
lithic core structure and low dissipation (Q - 2000) are as¬ 
sumed. The reason for this discrepancy is the difference in as¬ 
sumption of the asteroid internal structure. Moreover, its value 
could vary depending on inhomogeneity of the internal struc¬ 
ture. To consider the diffusion of seismic energy density, the 
highly inhomogeneous structure must be assumed. Here, we 
temporarily employ the constant value as standard ones mainly 
on the basis of Richardson et al. (2004); Richardson Jr et al. 
(2005). However, the details of parameter dependences will be 
discussed in section 4.2.3. The standard parameter values are 
summarized in Table 1. 

Then, we compute the duration of global seismic shaking Ty 
which is defined by a period of F > 1 in the entire asteroid. The 
specific form of Ty is written as. 


Tv(Di, A) = Tatten(T>i, A) “ Tdiff(A), (4) 


where Tatten(7Ji, TJa) is the period from the moment of impact to 
the time at which F decays to be 1 due to the inelastic dissipa¬ 
tion. And, Tdiff(Z)a) represents the time required to diffuse the 
seismic energy density to the entire asteroid. 



Figure 3: (a) Temporal decay of the maximum dimensionless acceleration F 
after the impact. The solid line shows an example of the case of Z); =4 m. 
Da = 400 m, and vi = 5.3, km To compute F(0, we use the standard 
parameter values: Q = 2000, rj = 10““*, f = 20 Hz, p[ = 1900 kg m“^, and 
pa = 2500 kg m“^. The horizontal dashed line con'esponds to F = 1. Vertical 
broken lines show and Tatten. and a two-way arrow represents Ty defined 
by Eq. (4). (b) Temporal decay of the convective velocity Vc. The solid cuiwe 
is computed by Eqs. (3) and (7) with the standard parameter values {Table 1). 
The hatched area corresponds to the migration length per impact / = 0.13 m. 
Its specific form is written in Eq. (9). 
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Table i: Standard parameter values for the global seismic shaking of small bodies. Mean impact velocity: v\, seismic quality factor: Q, impact seismic efficiency: 
7 ], dominant seismic frequency: /, density of the tai'get asteroid: pa, density of impactor: pj, and diffusion constant of seismic energy density: Kg. 


Vi (kms *) 

Q 

7 

/(Hz) 

P, (kg m 3) 

Pa (kg m ’) 

K, (m^ s ') 

5.3 

2000 

lo-'* 

20 

1900 

2500 

2.5 x 10^ 


The form of Tatten can be computed from Eq. (3) as, 


Tatten(A,T>a) = 

2nf 


3/vi PiDl 


plDl 


(5) 


The value of Tatten depends on both D\ and D^. In Fig. 4, 
Tatten(T*a) with Di - 400 m (pink solid line), Dt = 40 m (green 
solid line), and Dj = 4 m (purple solid line) are shown. All 
other parameter values are fixed to the standard ones in these 
lines. 

On the other hand, Tditf is the diffusion timescale of the seis¬ 
mic energy density (Richardson et al. (2004) and Richardson Jr 
et al. (2005)). The computed form of Tditf is written as, 

Dl 

7difF(T>a) - -7r~2- (6) 



Da [m] 


In Fig. 4, Tdiff(Z)a) is represented by the dashed curve. Here, 
we employ the typical value for the diffusion coefficient - 
2.5 X 10^ m^ s“^ (Richardson et al., 2004; Richardson Jr et al., 
2005). The global seismic shaking is achieved in the area above 
the dashed curve (Tdifr(Z)a)) and below the Tatten(JJi, JJa)- For 
instance, the two-way arrow in Fig. 4 shows Tv for the collision 
of Di = 4 m to Da = 400 m. The corresponding durations Tatten, 
Tditf and Tv are also shown in Fig. 3. 

If Tatten is Smaller than Tditf, the global seismic shaking is 
impossible. In such a situation, only the local seismic shaking 
is achieved in the vicinity of the impact point. While the global 
seismic shaking is attained at f = Tditf, the local seismic shaking 
occurs even in the early stage of the seismic shaking, t < Tditf. 
This is particularly the case for the large target asteroids of Da > 
5 X 10^ m. For the small target asteroids of Da < 5 x 10^ m, on 
the other hand, Tditf is negligibly short as shown in Fig. 4. For 
example, in the case of Itokawa-sized asteroid, the estimated 
diffusion timescale Tditf (Da = 400 m) = 6.6 x 10“^ s is much 
shorter than Tatten(Di = 4 m, Da = 400 m) = 180 s. In this 
regime, the global seismic shaking is quickly achieved, and 
can be simply approximated by Tatten- According to Fig. 4, the 
asteroids in the range of 5 x 10^ m < Da < 2 x 10^ m can be 
vibrated globally only for relatively a short time. Moreover, the 
asteroids in the range of Da > 2 x 10^ m cannot be vibrated 
globally because Tditf becomes longer than Tatten by the largest 
impactor which disrupts the target. (The chain line in Fig. 4; 
the definition of disruption limit is explained later in the next 
subsection.) Even for such large asteroids, the local seismic 
shaking could be induced at the neighborhood of the impact 
point. 

2.1.3. Regolith convection 

Substituting Eq. (3) into the velocity scaling of the granular 
convection (Eq. (8) in Yamada and Katsuragi (2014)), we can 


Figure 4: The relationships among various characteristic times for MBAs. The 
solid lines indicate the attenuating duration Tatten by the impactor with diameter 
Di = 400 m (pink), 40 m (green), and 4 m (purple) vs. the tai'get asteroid 
diameter D-^. The parameter values used are identical to those in Fig. 3. The 
dashed curve shows Tdiff with = 0.25 km“ s“^ (Richardson et al., 2004; 
Richardson Jr et al., 2005; Miyamoto et al., 2007). The chain curve shows Tatten 
by the largest impactor which disrupts itself (Jutzi et al., 2010). The global 
seismic shaking is attained in the hatched area. A two-way arrow shows the 
vibrated duration Ty = Tatten ” for the collision between the impactor of 
Di = 4m and target of D^ = 400 m (Itokawa size). 


estimate the velocity of regolith convection. Using the convec¬ 
tive velocity scaling and Tv, the migration length per impact / 
can be computed. Schematic image of / is represented by red 
arrows in Fig. 1. 

The migration length per impact / can be estimated as fol¬ 
lows. The granular convective velocity Vc has been character¬ 
ized as a function of T in the previous study (Yamada and Kat¬ 
suragi, 2014). The form of vdt. A, is written as, 


v^it, Di, Da) = Co 


/ r(f,Di,Da) y 

I 2nf I 



(7) 


where d is the diameter of grains and A is the convective roll 
scale. The roll size A is defined by the geometric mean of the 
horizontal and vertical dimensions A - 'JrH, where R and H 
are the radius and height of the vibrated granular bed. Yamada 
and Katsuragi (2014) determined the parameter values by the 
fitting to the experimental data as Co = 3.6 x 10“^, a = 0.47, 
and p - 0.82. We can obtain the following form of /(Di, Da) by 
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integrating Eq. (7) from Tdiff to Tatten: 


KA,£>a) 


^^attcn 

"Tdiff 


vc (r(f, A, A)) dt 


2Q'(27r/)i+2a 



f3fvi 

1 

2a 

X 

c 




( 8 ) 




Fig 3(b) shows Vc(f) calculated from Eqs. (3) and (7). In 
Fig. 3(b), the hatched area corresponds to 1(A = 4m, Da = 
400m) = 0.13 m. The values of Co, a, and p were obtained on 
the basis of spherical glass-beads experiments performed on the 
ground. Although it might not be appropriate to use these val¬ 
ues for regolith migration, here we employ these values for the 
first step approximation. Furthermore, parameter sensitivities 
in the estimate of resurface timescale will be discussed later in 
section 4.2 

To estimate the specific values of Vc and / using Eqs. (7) and 
(9), we have to assume the values of d and A in addition to Vi, Q, 
j], /, Pi, and Pa. We use 1 cm for the typical value of d based on 
the actual observation of Itokawa (Fujiwara et al., 2006; Yano 
et al., 2006). Although d obeys a certain distribution form in 
general, we use the typical constant value to investigate the rep¬ 
resentative behavior of the asteroid. For the convective roll size 
A, we assume that A fulfills A/c/ = 100. This assumption is 
based on the experimental observation of granular convection. 
In Yamada and Katsuragi (2014), a transition from one pair of 
convective rolls to two pairs of convective rolls was observed 
when the convective (horizontal) system size is approximately 
Ajd ^ 100. Aoki et al. (1996) also confirmed the splitting of 
convective rolls Ajd ^ 90. Furthermore, the convective roll 
most likely has a vertical intrinsic length scale as well. In the 
deeper part of the vibrated granular bed, the convective veloc¬ 
ity becomes almost zero. This quiescent part is called frozen 
zone (Taguchi, 1992; Ehrichs et al., 1995; Knight et al., 1996; 
Yamada and Katsuragi, 2014). From these experimental obser¬ 
vations, we assume that the granular convective roll size satis¬ 
fies Ajd - 100 both in horizontal and vertical directions (i.e., 
aspect ratio is 1) in the current model. 

Here, we employ additional restrictions to the range of Di 
that is relevant to compute the resurfacing timescale. Both the 
upper and lower limits for Di are introduced by considering the 
disruption and well-ordered convective motion by intermittent 
events, respectively. 

Too large impactors would disrupt the target asteroid. This 
upper limit value is denoted by Di ^ax- can be computed 

from the threshold of catastrophic disruption which corre¬ 
sponds to the critical amount of energy per unit mass deliv¬ 
ered by the impactor. The catastrophic disruption means that 
the target asteroid is pulverized and dispersed until the mass of 
the largest remnant is smaller than a half of the original mass. 
Q*^ for broad-size-ranging basalt targets was self-consistently 
characterized by Benz and Asphaug (1999) using numerical 
method. Jutzi et al. (2010) used an improved numerical code, 
so that Q*j for pumice targets, which are not fractured but mono¬ 


lithic with microporosity, can be computed. In addition, it 
should be noted that all Q*^ values introduced here are com¬ 
puted with continuum target modeling. In general, it is not easy 
to estimate the realistic Q*^ for rubble-pile targets. Since we 
focus on the porous rubble-pile or very fractured asteroids in 
this study, the porous-target-based Q*^ reported by Jutzi et al. 
(2010) is employed. The obtained Di,niax(A) is plotted by the 
chain curve in Fig. 5. Tatten by the disruptive impactor can also 
be computed from D^max- The chain curve in Fig. 4 indicates 
this disruption limit Tatten(T^i,max, Df) above which the target as¬ 
teroid is pulverized. 

On the other hand, too small impactor is not able to trigger 
the convective motion. In other words, I must be large enough 
to produce the well-ordered (reproducible) convective motion 
of the roll size A. Thus, we assume that the minimum limit 
of migration length /min should be related to the convective roll 
size A. If /min is much shorter than A, it is probably difficult 
to follow the same convective structure by a series of intermit¬ 
tent seismic shakings. Therefore, the lower limit Di ^in can be 
written as, 

Ain — uA, (10) 

where n is a parameter satisfying n > 0. The obtained lower 
limit Di_min(7Ja) in the case of n = 0.1 is presented by the solid 
curve in Fig. 5. Consequently, we hnd that there is a hnite area 
in which the global convective motion can be induced (hatched 
area in Fig. 5). Hereafter, we use n - 0.1 unless otherwise 
noted. Since we assume Aid - 100, n - 0.1 corresponds to 
Ain = 10(7. Although this n value is speculatively assumed, we 
evaluate its effect to the estimate of timescale by means of scal¬ 
ing method, later in section 4.2.3. The gray line in Fig. 5 corre¬ 
sponds to Di,min(7Ja) Computed with the approximation Tdifr = 0. 

Namely, Fig. 5 shows the phase diagram of the outcome in¬ 
duced by the impact between Di and Da. The hatched area in 
Fig. 5 represents the global seismic shaking phase that results 
in the well-ordered convective motion. Above the chain line 
(Di max), the target asteroid is destroyed by the impact. Below 
the solid or gray lines (Di ^in), the well-ordered convective mo¬ 
tion cannot be induced. This limit corresponds to the bound¬ 
ary for the no-convection phase due to the too small impactor. 
The upper right part bounded by chain, solid, and gray lines 
indicates the phase in which only the local seismic shaking is 
achieved. For Itokawa-sized asteroid (Da = 400 m), Di max and 
TJi.min are 10 and 3.4 m, respectively. 

In Fig. 6(a), the ideal convective velocities at the impact mo¬ 
ment (t = 0) are shown by three colored lines (Di = 400 m 
(pink), 40 m (green), and 4 m (purple)). The hatched areas 
correspond to the global seismic shaking regime. As shown in 
Fig. 6(a), the typical velocity for small asteroids is in the order 
of 10“^ m/s. This convective velocity is considerably smaller 

thantheescape velocity of the target asteroid Vesc = -^|7rGpaA 

as long as Da is greater than 10^ m. Thus, regolith grains would 
not disperse to space by regolith convection while a few grains 
may be released from the surface to space due to the spalla¬ 
tion. In addition to the spallation at the very surface, surface 
regolith layer might experience the long-time lofting (Garcia 
et ah, 2015). In such a situation, the convective migration might 
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(a) 



Da M 


Figure 5: A range of impactor size which can generate the regolith convec¬ 
tion of MBAs. The chain line represents the upper limit Z)i,max at which the 
disruption of the target asteroid occurs (Jutzi et al., 2010). The solid line repre¬ 
sents the lower limit Z)i,niin which is defined by the minimum migration length 
(n = 0.1 in Eq. (10)). The gray line shows Z)i,min computed with the approxi¬ 
mation Tdiff = 0. 


be significantly reduced. However, the quantitative modeling 
of this effect is still difficult since the details of the lofting (e.g., 
how thick it is, how long it lasts depending on vibration condi¬ 
tions) have not yet been revealed well. Therefore, we neglect 
the lofting effect in the current model. However, its quantitative 
modeling is one of the most important future steps to improve 
the model. 

The corresponding convective migration length per impact is 
shown in Fig. 6(b). The color and line codes used in Fig. 6(b) 
are identical to Fig. 6(a). The estimated maximum value of I 
is about 1 m, which coincides with the assumed convective roll 
size. Thus, the convective resurfacing can almost be completed 
by a single impact event if the impactor size is sufficiently large 
Di - £>i,max- However, such a large scale impact rarely occurs. 
As discussed later in section 3.2, the accumulation of small mi¬ 
grations is actually more effective than the single large impact. 
To compare the contributions by small and large impactors, the 
impact frequency Np has to be multiplied to 1. 


2.2. Resurfacing timescale 


Finally, by unifying all the steps mentioned above, we can 
obtain the resurfacing timescale as a function of the target as¬ 
teroid size, T(Df). Dividing the convective roll size A by the 
mean migration length per year L{Df), T(Df) can be calculated 


as. 


T{D,) 


A 

L(D,y 


( 11 ) 


Here, L(Df) can be computed by integrating the product of 
l{Di, Da) and Np{Di, Df). This procedure is described as. 


L{D,)^ 2 l(Di,Da)Vp(D.,Da). (12) 

A-F>i,min 



Da [m] 


Figure 6: (a) Regolith-convection velocity Vc on MBAs at the impact moment 
{t = 0). The chain curve represents the maximum convective velocity Vc,max 
which can be generated by the catastrophic impact {D\ = Z)i,max)- The solid 
curve represents the minimum convective velocity Vc,min t>y Di = Z)i,miiT The 
colored solid lines indicate Vc by the impactor with diameter Di = 400 m (pink), 
40 m (green), and 4 m (purple) vs. the tai'get asteroid of diameter D^. The red 
solid line represents surface escape velocity Vesc on the target asteroids. The 
horizontal gray line shows Vc,min computed with Z)i,min using the approxima¬ 
tion Tdiff = 0. (b) Migration length per impact / of MBAs. The chain cuiwe 
represents the maximum of / by the maximum impactor Di = A, max- The hor¬ 
izontal solid line represents the minimum limit of migration length The 
colored solid lines indicate I by various size impactors. The gray chain curve 
shows /max Computed with the approximation r^^ff = 0. The parameter values 
used in (a) and (b) are identical to those in Fig. 3. 
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The range of numerical integration in Eq. (12) spreads from 
to Z)i,niax- To compute the Specific value of L{D^) by 
Eq. (12), Eqs. (2), (9), and the data shown in Eig. 2 are used. 


3. Result and analysis 

3.1. Estimation of the resurfacing timescale 

The numerically computed T{Df) on the basis of OBG pop¬ 
ulation data for MBAs is shown by red solid circles in Eig. 7(a). 
Similarly, T{Df) for NEAs is shown by blue solid triangles in 
Eig. 7(b). Eor small asteroids of Da < 5 x 10^ m, the resurfacing 
timescale approximately obeys power law since the exponen¬ 
tial factor in Eq. (9) is negligible (almost unity due to the very 
short Tdiff). On the other hand, T(Df) tends to diverge in the 
large DJg> 5 x 10^ m) regime because the exponential factor in 
Eq. (9) dominates the behavior due to the large Tdiff. 

The mean collisional lifetime is also shown by black squares 
in Eigs. 7(a) and (b) (Jutzi et al., 2010; O’Brien and Green¬ 
berg, 2005). We can confirm that T{Df) necessary to resur¬ 
face the asteroid is shorter than the mean collisional lifetime 
for both MBAs and NEAs, as long as the order of asteroid sizes 
is smaller than 10^ m. Specifically, the resurfacing timescale 
of Itokawa-sized asteroid T (Da = 400 m) = 40 Myr is shorter 
than its collisional lifetime 1.7 x 10^ Myr. These results are in¬ 
sensitive to the difference between OBG and BAL populations. 

Even for the large asteroid of Da > 10“^ m in which the global 
seismic shaking is difficult, the resurfacing timescale might re¬ 
main short if we assume the local convection can contribute to 
the resurfacing process. The estimated T{Df) including local 
convection is shown by open red circles (MBA) or open blue 
triangles (NEA) in Eig. 7. These results (open circles and tri¬ 
angles in Eig. 7) have been reported in our previous conference 
abstract (Yamada et al., 2015). Actually, both estimates (with 
and without local convection cases) agree well in the small Da 
regime since Tdiff is sufficiently small to neglect the exponential 
factor in Eq. (9). 

However, we should notice that the obtained T (Da) can vary 
depending on parameters that are actually very uncertain. Eor 
example, the dominant frequency / might be larger than the 
assumed standard one particularly in small asteroids (Garcia 
et al., 2015). To verify the parameter dependences, the approxi¬ 
mated form of ^(Da) will be useful. Therefore, we compute the 
approximated analytic form of T (Da) in the next subsection. 

3.2. Scaling of the resurfacing timescale 

To obtain the approximated scaling form of T (Da), we begin 
with the power-law approximation of M,cum(^ D') for MBAs 
and NEAs as, 

M,cum(> D[) = C,D'r\ (13) 

where the value of Ci is a fitting parameter having dimensions 
of (lengthy and y is the characteristic exponent which corre¬ 
sponds to the slope of the solid line in Eig. 2. Cf for MBAs 
and NEAs are obtained as 3.2 x 10'^ m^-* and 5.2 x 10" m^'^. 



Da [m] 

Figure 7: (a) Various timescales for MBAs. Red solid circles represent the 
timescale necessary to resurface the asteroid by regolith convection, T{Da), 
which is numerically computed by Eq. (11). Red open circles represent T{Da) 
by considering local convective motion as well. The red solid line shows the 
approximated T{Da) computed by Eq. (20). Black solid squares represent the 
mean collisional lifetime for MBAs (Jutzi et al., 2010; O’Brien and Green¬ 
berg, 2005). Population data for MBAs (red circles in Fig. 2) are used for the 
computation. An error bai' at the Itokawa size indicates the cosmic-ray expo¬ 
sure (CRE) age of Itokawa’s returned particle (Nagao et al., 2011; Meier et al., 
2014). T{Da) are computed by using the standard parameter values, (b) Various 
timescales for NEAs. Blue solid triangles represent T{Da) for NEAs computed 
by Eq. (11). Blue open triangles represent T{Da) by considering local con¬ 
vective motion as well. The blue solid line shows the approximated T{Da) 
computed by Eq. (21). Population data for NEAs (blue triangles in Fig. 2) are 
used. T{Da) are computed by using the standard parameter values. The blue 
band corresponds to orbital lifetime for NEAs (Gladman et al., 2000) 
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respectively, y for MBAs and NBAs are 2.8 and 3.2, respec¬ 
tively. Equations (2), (12), and (13) enable us to analytically 
approximate L{D^) as. 


B(£>a) 



l(Di,D,)Np{Di,D,) 


dDi 




DilnB 

A A'" 




I £(M,cum(>A)\ dDi 
\ c/A/AlnB JDilnB’ 


(14) 


(15) 


where the minus sign in front of differential term of Ai_cum 
is for the negative slope of M,cum(^ D') (i.e., y > 0), and 
B = 10°' ^ 1.3 is the constant ratio between D', and 
The factor InB comes from the logarithmic binning effect, 
dloggDi - dDilDilnB (O’Brien and Greenberg, 2005). Sub¬ 
stituting Eqs. (9) and (13) into Eq. (15), the integration can be 
performed analytically. 

Eq. (15) yields twelve terms that depend on A,max or Z)i,min- 
By comparing the orders of these terms, a principal term domi¬ 
nating the behavior of L{D^ can be identified. Using the dom¬ 
inant term, L{D^) is approximated as. 


AA) 


Q, ~-2a+\ „ ^ 

XPiPa ~PiCl 


xD/“''^D^“T^exp 

i.min ^ 


Aanf \ 

-^TdifrJ, 


where C* is a dimensionless constant. 


(16) 


rCo(|;r)“^"32“ 

4(3a-r)2Q'(27r)2“+^' 

In Eq. (16), the dominant term depends on Di min. This implies 
that the small impactors play an essential role for the regolith 
resurfacing process. Since Z)i,min is determined by Ain through 
Eq. (9), a little more computation is necessary to obtain the final 
form for the analytic approximation of L{D^). Then, T{D^) can 
be readily calculated by Eq. (11). 

Erom the above approximation and a little more algebra, we 
finally obtain an approximated scaling form of ^(Da) expressed 
by the product of simple power-law forms and an exponential 
form. The obtained form is written as. 


T(A) 


„ -Z J.JL ,2 + &_JL 

CttQ ^“rj 6a 

-zL y By T 7 -Z 2Z--L 




Ay 

3 

a 


exp 


/ 4y7r/Tdiff \ 

I 30 j’ 


(18) 


where C„ is a dimensionless constant. 


C„ 


A (3a - y) 
7 


CQ{jn) 


2„)“+2 32 a 


2a (In) 


2 ( 7+1 



(19) 


Substituting a - 0.47, yS = 0.82 (experimentally obtained 
values), and y = 2.8, the specific form of Eq. (18) for MBAs is 


obtained as. 


r(A) = 


gl.97^0.93^,l.S5^0.87p.(50.06Cj 


£>a''exp 


3.77r/Tdiff \ 

Q /’ 

( 20 ) 


where Ci = 3.2 x 10'^ m^ * is used for MBAs. The solid red 
curve in Eig. 7(a) represents Eq. (20). Using the identical a and 
P values, and y = 3.2, Eq. (18) for NBAs is also computed as. 


T{D,) 


C„n‘-23A°-4V2-23^i.76p0.98 /4.2;r/Td.ff \ 

g2.23^1.05^,2.10^1.05p.(50.07Ci ^ Q j’ 

( 21 ) 


where Ci = 5.2 x 10" m^-^ is used for NBAs. The solid blue 
curve in Eig. 7(b) represents Eq. (21). We can confirm that both 
approximated forms of ^(Da) agree very well with numerically 
computed ones as long as we neglect the local seismic shaking 
effect. Moreover, Eqs. (20) and (21) allow us to evaluate the 
various parameter dependences of T (Da) for MBAs and NBAs, 
respectively. 


4. Discussion 

4.1. Comparison with other timescales 

We think that the asteroids have been resurfaced during their 
stay in the main belt. Since the orbital lifetime for NBAs is 
ranging from 1 Myr to 10 Myr (Gladman et al., 2000), ^(Da) 
for NBAs in Eig. 7(b) is greater than the orbital lifetime. This 
means that it is difficult to resurface the asteroids within their 
orbital lifetime for NBAs. Thus, we basically consider that the 
resurfacing is important for MBAs, and discuss hereafter on 
MBAs values although both of Eros and Itokawa currently be¬ 
long to NBAs. 

In Eig. 7, the mean collisional lifetime (for both MBAs and 
NBAs) is longer than T{Dp. However, the mean collisional 
lifetime depends on the model of Q^. As mentioned in Sec¬ 
tion 2.1.3, the rubble-pile asteroid is more fragile than the 
monolith asteroid. Thus, for rubble-pile asteroid might be 
much smaller than the estimate by Jutzi et al. (2010) (e.g., 
Durda et al. (1998)). Such a small actually results in the 
short collisional lifetime. Namely, the uncertainty of also 
affects the competition between the mean collisional lifetime 
and T(Dp. Note that, however, Q*^ does not directly affect the 
estimate of TfDa) because it is rather dominated by Di min. 

Eor asteroid Itokawa, Michel et al. (2009) suggested that the 
timescale necessary to reproduce the current crater record (Hi- 
rata et al., 2009) ranges from 75 Myr to 1 Gyr. To estimate the 
timescale, they used the impactor population model (O’Brien 
and Greenberg, 2005; Bottke Jr et al., 2005), seismic shaking 
model (Richardson et al., 2004; Richardson Jr et al., 2005), and 
crater scaling law (Holsapple, 1993; Nolan et al., 1996). The 
estimated timescale corresponds to the surface age computed 
from the number density of craters. They found that the surface 
age of Itokawa is equivalent to or larger than the corresponding 
collisional lifetime in the main belt. In addition, the resurfacing 
timescale of Itokawa is much shorter than this crater-based sur¬ 
face age. This suggests that the regolith convection may keep 
the surface of Itokawa fresh by the resurfacing during the crater 


9 
















erasure process. Moreover, the characteristic surface structure 
such as impact craters might be erased by the regolith convec¬ 
tion. Then, the relatively small craters of diameter smaller than 
10 m might be erased on the surface of small asteroids because 
the assumed convective roll size is in the order of 10*’ m. In¬ 
deed, the number density of craters smaller than 10 m shows an 
unexplained decrease on Itokawa’s surface (Hirata et al., 2009; 
Michel et al., 2009). 

The value of resurfacing timescale on Itokawa 40 Myr is 
longer than its CRE age (1.5-8 Myr) (Nagao et al., 2011; Meier 
et al., 2014). If we assume that the regolith grains have been ir¬ 
radiated with cosmic rays while they are on the surface (upper 
part) of Itokawa in the course of convective motion, the resur¬ 
facing timescale should be equal to or longer than the CRE age. 
Thus, this relatively young CRE age is consistent with the cur¬ 
rent regolith convection model. 

Eor asteroid Eros, Richardson et al. (2004) and O’Brien et al. 
(2006) estimated the surface age to reproduce the observed 
crater record reported by Chapman et al. (2002). The esti¬ 
mated surface age for Eros is 400 + 200 Myr (Richardson et al., 
2004) or 120 Myr (O’Brien et al., 2006). Although the differ¬ 
ence between these two ages originates from the difference of 
the used crater scaling laws, both ages are much shorter than 
the mean collisional lifetime (8.9 x 10^ Myr) in MBA. These 
estimated surface ages are far below the resurfacing timescale 
T{D^ = 1.7 X 10"^ m) = 6.6 x 10^ Myr. And, this resurfacing 
timescale ( 6.6 x 10^ Myr) is much longer than the mean colli¬ 
sional lifetime as well. This means that it is difficult to resur¬ 
face the asteroid Eros by global regolith convection. However, 
by considering the effect of the local regolith convection in the 
vicinity of impact point (open circles in Eig. 7(a)), T{D^) for 
asteroid Eros (5.2 x 10^ Myr) would be close to these estimated 
surface ages. 

4.2. The evaluation of resurfacing feasibility 

As already mentioned in the last paragraph of Section 3.1, 
T{Df) involves a number of uncertain parameters. Thus, the ef¬ 
fect of parameter uncertainties to the feasibility of resurfacing 
by regolith convection should be discussed. To evaluate the un¬ 
certainty effect, we first reconfirm the range of each parameter 
uncertainty in Section 4.2.1. Then, we introduce a dimension¬ 
less coupling parameter to simply parametrize the possibility 
of regolith convection (Section 4.2.2). Einally, the feasibility of 
convective resurfacing is discussed in the parameter space and 
the dimensionless coupling parameter (Section 4.2.3). 

4.2.1. The ranges of parameter uncertainty 

Eqs. (20) and (21) allow us to read the relative parameter 
sensitivities in the form of T(Df) for MBAs and NEAs, respec¬ 
tively. The larger exponent value implies the stronger depen¬ 
dency in each factor of Eqs. (20) and (21). In this sense, Q, 
/, d, and Vi are relatively sensitive parameters, whereas p^, pi, 
n, A, and 77 are relatively insensitive parameters. Actually, the 
values of Q, /, n. A, and 77 are very uncertain. Ranges of these 
uncertain parameters are tabulated in Table 2. The values of d, 
v'i. Pa and Pi can somewhat be estimated by observation and/or 


appropriate modelings. Thus, the accurate estimate of Q and / 
is the most important factor for the reliable estimate of T{Df). 
For the quality factor Q, we consider that its order ranges from 
200 to 2000 by following Richardson et al. (2004) and Richard¬ 
son Jr et al. (2005). We consider that S-type asteroids should 
not be dissipative due to the lack of dissipative components. 
Thus, Q - 2000 would be more plausible. For the worse case 
estimate, we consider the case Q - 200 as well. The domi¬ 
nant seismic frequency computed by Richardson et al. (2004) 
and Richardson Jr et al. (2005) is / = 10-20 Hz. However, the 
recent numerical study revealed that / could be greater than 20 
Hz (perhaps over 100 Hz) particularly for small asteroids (Gar¬ 
cia et al., 2015). The value of / seems to strongly depend on 
the internal structure of target asteroids. Therefore, further sys¬ 
tematic investigation is necessary to estimate the reliable value 
of the dominant vibrational frequency. Here, we assume that / 
might vary from 10 to 200 Hz. The asteroids considered in this 
study is highly inhomogeneous with low dissipation. 


Table 2: Lower and upper limits of parameters used in the model. 



lower limit 

standard 

upper limit 

Q 

200 

2000 

2000 

/(Hz) 

10 

20 

200 

7 

10 -® 

10-4 

10 ' 

Aid 


100 


n 

0.001 

0.1 

0.1 


The value of 77 ranges from 10“^ to 10 ' (Richardson et al., 
2004; Richardson Jr et al., 2005). Because the uncertainty of 77 
spreads over five orders of magnitude, it may also significantly 
affect the estimate of T(Df). Yasui et al. (2015) have recently 
reported that 77 ranges from 5 x 10“^ to 5 x 10 Although 
this result is based only on the relatively low impact velocity 
experiment, the estimated value is consistent with the assumed 
standard value. 

In addition, A and n are particular parameters in the present 
model. For the convective roll size A, we assume A/t/ = 100 on 
the basis of experimental observations (Yamada and Katsuragi, 
2014; Aoki et al., 1996). However, the physical condition de¬ 
termining the intrinsic convective roll size has not yet been re¬ 
vealed. As for 77 , while we temporarily assume n - 0.1, the rea¬ 
son for employing this value is not firmly based. The basis of 
this speculation (77 = 0 . 1 ) is the preliminary observation of the 
tapped granular layer (likawa et al., 2015). According to the ex¬ 
perimental result, however, very small /min(- O.OOIA) may also 
be able to produce a well-ordered convective structure by the 
accumulation of intermittent tappings. Since T{Df) is propor¬ 
tional to 77 ° (Eq. (20)) or 77 ' (Eq. (21)), the smaller n results 
in the shorter T{Df). This tendency does not significantly af¬ 
fect the interpretation of the current result, viz., we employ the 
stricter condition for n in the current estimate. However, note 
that the experimental observation is still quite preliminary. The 
accurate determination of A and n is a very important future 
problem to be investigated. 

As discussed in Section 1, gravity dependence of granu¬ 
lar convection is still under debate. Although here we em¬ 
ploy the values of Cq, a, and f obtained by Yamada and Kat- 
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suragi (2014), these values are also uncertain and can affect the 
timescale estimate (Eqs. (18) and (19)). 


4.2.2. Dimensionless parameter for resurfacing difficulty 
Next, we can introduce a dimensionless coupling parameter 
for evaluating the feasibility of asteroid resurfacing by regolith 
convection. In the following calculation, we use MBAs value 
y = 2.8 as well as a = 0.47 and /3 - 0.82 to clearly observe the 
parameter dependences. Since the form for NBAs is basically 
identical to that for MBAs, the conversion to the form for NBAs 
is straightforward. By neglecting the exponential factor, a non- 
dimensionalized form of Bq. (20) can be described as. 


Y 


/ 




2.78\-1.02 


Gp.Dl 


0.44 



( 22 ) 

where F is a dimensionless parameter which consists of princi¬ 
pal uncertain parameters, Q, rj, n, /, d, and A. Then, the specific 
form of Y can be obtained as. 


Y = C„Q 


-1.97 -0.93„0.97 


/ 


ViCj 


0.36 


1.97 




0.36x0.35 


(23) 

This dimensionless parameter Y indicates the difficulty of re¬ 
golith convection. Approximated T{Df) shown in Big. 7(a) (the 
red curve) corresponds to the case of T ^ 10“^. Bor NBAs, 
the blue curve in Big. 7(b) corresponds to T ^ 10^^. Actually, 
T{Df) with T ^ 10-"^ for MBAs (or Y ^ 10^® for NBAs) be¬ 
comes comparable to the mean collisional lifetime for MBAs 
(or NBAs). This implies that it is difficult to resurface the aster¬ 
oid with - log 10 Y < 4 for MBAs (or - logio T < 6 for NBAs) 
because such an asteroid would be disrupted by the fatal col¬ 
lision before its resurfacing. If - logiQ Y is greater than 4 for 
MBAs (or 6 for NBAs), the regolith convection can be a rele¬ 
vant process for resurfacing the asteroids. Using the - logjQ Y 
value, we can quickly evaluate the plausibility of the resurfac¬ 
ing by regolith convection once the accurate parameter values 
are obtained by explorations, experiments, numerical model¬ 
ings, and so on. 


Possible 

Impossible 

d =10 ^ m 



f = 20Hz\,.^^ 

d=10 m 

f = 200 Hz! 

• Possible 




'4 






10 



f[Hz] 


4.2.3. Scaled parameters to resurfacing 

To systematically evaluate the possibility of asteroid resur¬ 
facing by regolith convection, we use Bq. (23) and the crite¬ 
rion - logio T > 4 in parameter spaces. As discussed in sec¬ 
tion 4.2.1, Q, /, and d are very sensitive parameters to estimate 
the resurfacing timescale. In addition, 77 has an extremely wide 
range of uncertainty. Thus, we focus on these parameter depen¬ 
dences in the following. In Big. 8 (a), the resurfacing threshold 

- logio T = 4 is drawn in the {d vs. /) parameter space. The 
threshold levels with Q - 2000 and 200 are then drawn by blue 
and green lines, respectively. Other parameter values are fixed 
to the standard ones, i.e., 77 = 10 is used in Big. 8 (a). In the 
region below the threshold level, the resurfacing can take place 
within the mean collisional lifetime. 

At the standard parameter point (solid circle in Big. 8 (a)), 

- logio ^ sufficiently larger than 4 because the threshold 
level is much above this point (blue line in Big. 8 (a)). The 


Figure 8: Resurfacing diagram in {d vs. /) parameter space. The square re¬ 
gion corresponds to the range of uncertainty of parameters (10~^ < d < 10“^ 
m (Fujiwara et al., 2006; Yano et al., 2006) and 10 < / < 200 Hz (Richard¬ 
son et al., 2004; Richardson Jr et al., 2005; Garcia et al., 2015)). The black 
solid circles represent the standai'd point (d = 10“^, / = 20 Hz). The colored 
solid lines represent the resurfacing threshold level (- logjQ Y = 4) at which 
the resurfacing timescale becomes comparable to the mean collisional lifetime. 
These lines can be computed from Eq. (23). The hatched area con'esponds to 
the resurfacing-possible regime (-logjQ F > 4) with respective parameter, (a) 
Q and (b) t]. Other parameters are fixed to the standard ones. 
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central square region in Fig. 8 (a) (10^^ < d < 10“^ mm x 
10 < / < 200 Hz) corresponds to the parameter uncertainty 
range considered in this study. We find that most of the square 
region is covered by the hatched area (regime below the thresh¬ 
old level with Q — 2000). This suggests that the resurfacing 
by regolith convection can occur in various S-type of asteroids. 
Even if the dominant frequency of seismic shaking is quite high 
(200 Hz) as reported by Garcia et al. (2015), small regolith 
grains {d ^ 10 “^ m) can be resurfaced within the mean colli- 
sional lifetime. However, there is a small regime at which the 
resurfacing timescale becomes longer than the mean collisional 
lifetime in Fig. 8 (a). The white triangular region at the top right 
corner in the uncertain square corresponds to this hardly resur¬ 
faced regime. 

Although the 77 does not have a strong sensitivity in Eq. (20), 
its value has an extremely wide range of uncertainty. To evalu¬ 
ate the effect of uncertainty, we draw another threshold levels in 
the identical parameter-space diagram in Fig. 8 (b). In Fig. 8 (b), 
the threshold levels (-logT = 4) are shown with rj - 10^', 
10 and 10“^. Other parameter values are fixed to the stan¬ 
dard ones, i.e., Q - 2000 is used in Fig. 8 (b). If 77 is sufficiently 
large, the resurfacing can easily occur within the mean colli¬ 
sional lifetime even if the seismic frequency is high (200 Hz). 
However, if the energy transfer is not so effective (77 = 10“^), 
the resurfacing timescale becomes longer than the mean colli¬ 
sional lifetime. 

As discussed above, it is difficult to conclude the possibil¬ 
ity of asteroid resurfacing only from the currently available in¬ 
formation. To accurately estimate the possibility, we have to 
restrict the parameter values in much smaller regime. In this 
study, we provide a simple framework by which the possibility 
of asteroid resurfacing by regolith convection can be evaluated 
once these parameter values were well determined. A stochas¬ 
tic modeling such as Monte Carlo simulation might be help¬ 
ful to evaluate the possibility. The stochastic modeling could 
also be useful to discuss the random accumulation of the re¬ 
golith migration composing a convective motion. These types 
of stochastic modelings are important future problems to refine 
the estimate reliability. 

4.3. Effects not considered in the model 

In the model, convective motion of regolith grains consists 
of the accumulation of intermittent migrations. However, the 
current form of granular convective velocity is obtained on the 
basis of steady vibrated granular experiment (Yamada and Kat- 
suragi, 2014). The granular convective velocity scaling under 
the truly intermittent attenuating vibration should be revealed 
to improve the model. The values of Co, a, and (3 might be 
affected by such experimental conditions. 

There are some other issues that might significantly affect the 
estimate of T{D^). For instance, the effect of grains cohesion 
might be essential in the regolith convection under the micro¬ 
gravity condition (Scheeres et al., 2010), although we neglect 
its effect in the current model. The effects of polydispersity 
and irregular shape of regolith grains are also omitted in this 
study. Some other factors that should be considered regarding 
experimental conditions such as container wall and interstitial 


air were also discussed in the previous paper (Yamada and Kat- 
suragi, 2014). In Garcia et al. (2015) or Scheeres et al. (2010), 
lofting of the regolith layer by impact-induced vibration was 
studied. During the lofting, regolith grains cannot make a con¬ 
tact network. Then, the regolith convection might hardly occur 
in very shallow part of the regolith layer. By the lofting ef¬ 
fect, the convective regime (hatched areas in Figs. 4, 5, and 6 ) 
could become smaller than the current estimate. More accu¬ 
rate and quantitative modeling of the lofting is a possible future 
work. The current model only discusses the considerably aver¬ 
aged picture of the global regolith convection. 

5. Conclusions 

We developed a simple model of regolith convection which 
enables us to estimate the resurfacing timescale caused by re¬ 
golith convection. The model consists of three steps: impact, 
vibration, and the accumulation of intermittent convections. Al¬ 
though we assume the extremely averaged picture of the as¬ 
teroid covered with regolith, we can estimate the resurfacing 
timescale of the various-sized S-type asteroids of diameter Da, 
using the model. To estimate the resurfacing timescale, the con¬ 
stant standard values are used although some of them depend 
on asteroid properties. In addition, we present a dimensionless 
parameter Y that indicates the difficulty of resurfacing process. 
Using Y, the rough estimate of the possibility of resurfacing can 
be performed. Because a number of uncertain parameters are 
used in the present model, the approximate scaling form of T 
and Y are helpful to evaluate the relative importance of the pa¬ 
rameters. As a result, we find that it is possible to resurface the 
asteroids by regolith convection within their mean collisional 
lifetime if we employ reasonable parameter value ranges. Thus, 
the regolith convection should be taken into account as one of 
the plausible resurfacing processes in the history of asteroid sur¬ 
face evolution. The knowledges obtained by further observa¬ 
tions and experiments will refine the regolith convection model. 
In particular, the observation by Hayabusa2 and OSIRIS-REx or 
their returned samples will bring us various helpful information 
to estimate much more reliable resurfacing timescale. 
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